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Inverted spectroscopy and interferometry for quantum-state reconstruction of systems 

with SU(2) symmetry 
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Department of Physics, Technion — Israel Institute of Technology, Haifa 32000, Israel 

We consider how the conventional spectroscopic and interferometric schemes can be rearranged 
to serve for reconstructing quantum states of physical systems possessing SU(2) symmetry. The 
discussed systems include a collection of two-level atoms, a two-mode quantized radiation field with 
a fixed total number of photons, and a single laser-cooled ion in a two-dimensional harmonic trap 
with a fixed total number of vibrational quanta. In the proposed rearrangement, the standard 
spectroscopic and interferometric experiments are inverted. Usually one measures an unknown 
frequency or phase shift using a system prepared in a known quantum state. Our aim is just 
the inverse one, i.e., to use a well-calibrated apparatus with known transformation parameters to 
measure unknown quantum states. 
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I. INTRODUCTION 



The last few years were marked by an outburst of re- 
search devoted to the problem of reconstruction of quan- 
tum states for various physical systems (see, e.g., Ref. ||] 
for an extensive list of the literature on the subject). The 
problem, as stated already in the fifties by Fano [g[ and 
Pauli 1^ , is to determine the density matrix p from infor- 
mation obtained by a set of measurements performed on 
an ensemble of identically prepared systems. Significant 
theoretical and experimental progress has been achieved 
during the last decade in the reconstruction of quantum 
states of the light field [Q. Also, numerous works were 
devoted to reconstruction methods for other physical sys- 
tems. Most recently, a general theory of quantum-state 
reconstruction for physical systems with Lie-group sym- 
metries was developed 

In the present work we consider state-reconstruction 
methods for some quantum systems possessing SU(2) 
symmetry. The principal procedure for the reconstruc- 
tion of spin states was recently presented by Agarwal . 
A similar approach was also proposed by Dodonov and 
Man'ko while the basic idea underlying this method 
goes back to the pioneering work by Royer 0|. In brief, 
one applies a phase-space displacement [specifically, a ro- 
tation in the SU(2) case] to the initial quantum state and 
then measures the probability to find the displaced sys- 
tem in a specific state (the so-called "quantum ruler" 
state). Repeating this procedure with identically pre- 
pared systems for many phase-space points [many rota- 
tion angles in the SU(2) case], one determines a function 
on the phase space (the so-called operational phase-space 
probability distribution |p[-|l0|). In particular, by mea- 
suring the population of the ground state, one obtains the 
so-called Q function. The information contained in the 
operational phase-space probability distribution is suf- 
ficient to completely reconstruct the unknown density 
matrix of the initial quantum state. A general group- 
theoretical description of this method and some exam- 



ples, including SU(2), are presented in Ref. 0]. 

The aim of the present paper is to study how the gen- 
eral state-reconstruction procedure outlined above can be 
implemented in practice for a number of specific physical 
systems with SU(2) symmetry. Three systems are consid- 
ered: a collection of two-level atoms, a two-mode quan- 
tized radiation field with a fixed total number of pho- 
tons, and a single laser-cooled ion in a two-dimensional 
harmonic trap with a fixed total number of vibrational 
quanta. We show that a simple rearrangement of conven- 
tional spectroscopic and interferometric schemes enables 
one to measure unknown quantum states of these sys- 
tems. 



II. RECONSTRUCTION OF QUANTUM STATES 
FOR SYSTEMS WITH SU(2) SYMMETRY 

We start with some basic properties of SU(2) which is 
the dynamical symmetry group for the angular momen- 
tum or spin and for many other systems (e.g., a collection 
of two-level atoms, the Stokes operators describing the 
polarization of the quantized light field, two light modes 
with a fixed total photon number, etc.). The su(2) simple 
Lie algebra consists of the three operators {Jx, Jy, Jz}i 



(1) 



The Casimir operator is a constant times the unit op- 
erator, — j{j + 1)1, for any unitary irreducible rep- 
resentation of the SU(2) group; so the representations 
are labeled by the single index j that takes the values 
j = 0,1/2,1,3/2,.... The representation Hilbert space 
Hj is spanned by the complete orthonormal basis fj,)} 
(where = j, j - 1, . . . , -j): 

J^b» + Jz\j,iJ.) ^ ii\j,ii). 

In the following we assume that the state Itp) of the sys- 
tem belongs to Hj (or, for mixed states, that the density 
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matrix p is an operator on Wj). Group elements can be 
parametrized using the Euler angles a, /9, 7: 



5 = g(a,/3,7)=e»-^^e'^'^«e^^-'^. 



(2) 



We will employ two very useful concepts: the phase 
space (which is the group coset space of maximum sym- 
metry) and the coherent states (each point of the phase 
space corresponds to a coherent state). For SU(2), the 
phase space is the unit sphere §^ = SU(2)/U(1), and each 
coherent state is characterized by a unit vector pl] , p^ 



n = (sin 6 cos sin 9 sin 0, cos i 



(3) 



Specifically, the coherent states |j;n) are given by the 
action of the group element 



5(n) = e-"t''''e-''^-'^ 



(4) 



on the highest-weight state 
|j;n) ^g{n)\j,j) 
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cosJ'+''(6'/2) 
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xsin^-''(0/2)e-^''*|j- m). (5) 



An important property of the coherent states is the res- 
olution of the identity: 



2i + i 

An 



dn|j;n)(j;n| 



(6) 



where dn = sin 9 d9 d(f>. 

A possible procedure for the quantum-state reconstruc- 
tion is as follows |]l|,^|J§]- First, the system, whose initial 
state is described by the density matrix p, is displaced in 
the phase space: 



p(n) =.g-i(n)p.9(n), n £ S^. 



(7) 



Then one measures the probability to find the displaced 
system in one of the states \ j, p) (e.g., in the highest state 
|j, j)). This probability 



P/^(n) = (j,M|p(n)|j» 



(8) 



(which is sometimes called the operational phase-space 
probability distribution) can be formally considered as 
the expectation value 



p^(n) =Tr[pr^(n)] 
of the so-called displaced projector 

rp(n) = g(n)|j»(j>|g"i(n). 



(9) 



(10) 



Repeating this procedure (with a large number of identi- 
cally prepared systems) for a large number of phase-space 
points n, one can determine the function p^(n). 

Knowledge of the function p^in) is sufhcient for the 
reconstruction of the initial density matrix p. We can 



use the following expansion for the density matrix (such 
an expansion exists for any operator on fij)'. 

2j I 

P = J2J2 T^imDim, Tlim = Tr {pDlj. (11) 

Here, Dim are the so-called tensor operators (also known 
in the context of angular momentum as the Fano multi- 
pole operators fl^ ). 



21 + 1 
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where (ji, mi; j2, m2|j, m) are the Clebsch-Gordan coef- 
ficients. Now, one can reconstruct the density matrix by 
using the relation 



v/(2j + I)/ An 
{j,tJ-'J,0\j,p) Js- 



f dnp^{n)Y*^{n), (13) 



where Yimin) are the spherical harmonics. Other ways 
to deduce the density matrix from the measured proba- 
bilities (n) were also proposed ||l^ . 

Let us also consider the useful concept of phase-space 
quasiprobability distributions (QPDs). In the SU(2) 
case, one can introduce an s-parametrized family of the 

QPDs 



^(n; ^) = E E ^/^MMt^Mu). (14) 



/=0 m=-l 



For s = 0, we have the SU(2) equivalent of the Wigncr 
function, 



Win) 
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;^E E '^ln^YlM. ( 
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1=0 m=-l 



For s = 1, we obtain the SU(2) equivalent of the Glauber- 
Sudarshan function (also known as Berezin's contravari- 
ant symbol), P(n), whose defining property is 



P^^^ I dnP(n)|j;n)(j;n|. 



(16) 



The function which is probably the most important for 
the reconstruction problem is the SU(2) equivalent of the 
Husimi function (also known as Berezin's covariant sym- 
bol), 



Q(n) = (j;n|p|j;n). 



(17) 



obtained for s = — 1. As is seen from Eq. (g), the function 
(5(n) gives the probabihty to find the displaced system 
in the highest spin state \ 



Q(n) =Pj(n). 



(18) 
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Also, one can see that the probability p-j {9, </>) to find the 
displaced system in the lowest spin state |j, —j) is equal 
to Q{9 + TT,(f)). More generally, any one of the QPDs can 
be reconstructed using the relation ^ 



P(n;s) = ^^^dn'/^-,(n,n')p^(n'). 



(19) 



2j 



1=0 



0',M;',o|j,M) 



where Pi{x) are the Legendre polynomials. For s = — 1 
and n = j we recover the relation (nS) . 



III. GENERAL DESCRIPTION OF 
EXPERIMENTAL SCHEMES 

A. Spectroscopy and interferometry 



where the transformation operator is 



_ „MlJy „iV>Jz „i-62Jy 



(24) 



Now, the aim is to measure the value of ip which is pro- 
portional to the transition frequency in a spectroscopic 
experiment or to the optical path difference between the 
two arms of an interferometer. The information on ip is 
inferred from the measurement of the observable Jz at 
the output. The quantum uncertainty in the estimation 
of ip is 



Aif ■■ 



AJz. 



\d{Jzout)/dip\'' 



(25) 



where the expectation values are taken over the initial 
quantum state of the system. This state is assumed to 
be known, so one can estimate the value of ip and the 
corresponding uncertainty. 



Quantum transformations which constitute the basic 
operations in spectroscopic and interferometric measure- 
ments can be conveniently described as rotations in an 
abstract 3-dimensional space. In this description, the 
system is characterized by the vector J = {Jx, Jy, Jz)"^ , 
where the three operators Jx, Jy, and Jz satisfy the su(2) 
algebra (||). 

A spectroscopic or interferometric process is usually 
described in the Heisenberg picture as a unitary trans- 
formation 

Jout = U{di,d2,ip)JUH'&U^h,^) = \J{^1,^2,'P)J, (21) 

where U('!9i, 192, v) is a 3x3 transformation (rotation) ma- 
trix, and i?2, if are transformation parameters (rota- 
tion angles). A standard transformation consists of three 
steps: 

(i) rotation around the y axis by i?!, with the trans- 
formation matrix Rj,('!9i), 

(ii) rotation around the z axis by ip, with the transfor- 
mation matrix ((/?), 

(iii) rotation around the y axis by 1^2 , with the trans- 
formation matrix Rj,(i?2)- 

The overall transformation performed on J is 

U{9,cb) = Ry{d2)Rz{'p)Ry{^l)- (22) 

This transformation is slightly more general than those 
routinely made in spectroscopy and interferometry. The 
usual choice is 1^2 = —^1 — ±7''/2, so U = Rx{±ip), re- 
spectively, while if is the parameter to be estimated in 
the experiment. In the Schrodinger picture, the density 
matrix of the system transforms as 

Pout = U\^i,d2,ip)pU{^i,^2,^), (23) 



B. Reconstruction of the initial state 

In this paper we consider how to use the spectroscopic 
or interferometric arrangement for the inverse purpose, 
i.e., for the measurement of an unknown initial quantum 
state by means of a large number of transformations with 
known parameters. 

As discussed in Sec. II, the first part of the recon- 
struction procedure is the phase-space displacement of 
Eq. (0). With the phase space being the sphere, this 
displacement is just a rotation produced by the opera- 
tor g{n) of Eq. (^). Now, compare this rotation with 
the one made during a spectroscopic or interferometric 
experiment, as given by Eqs. ( p3| ) and (p4|). One can 
immediately conclude that the phase-space displacement 
needed for the SU(2) state-reconstruction procedure can 
be neatly implemented by means of the spectroscopic and 
interferometric techniques. One only needs to omit the 
first rotation (i.e., take Si =0), and recognize the two 
spherical angles as: 

9 = -^2, (t> = (26) 

After the rotation g{r\) is made, one should measure the 
probability p^(n) to find the displaced system in the state 
I j, /i) . Perhaps the most convenient choice is to measure 
the population of the lowest state |j, — j), which is usu- 
ally the ground state of the system (e.g., this state corre- 
sponds to the case where all the atoms are unexcited; 
in the atomic case such a measurement can be made 
by monitoring the resonant fluorescence for an auxiliary 
dipole transition). This procedure should be repeated for 
many phase-space points n with a large number of identi- 
cally prepared systems, thereby determining the function 
p^(n) (e.g., for = j or /i = —j). According to the for- 
malism presented in Sec. II, this information is sufficient 
to reconstruct the initial quantum state. 
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IV. COLLECTIONS OF TWO-LEVEL ATOMS 

In Ramsey spectroscopy Q one deals with a collec- 
tion of N two-level systems (usually atoms or ions) in- 
teracting with classical light fields. One can equivalently 
describe this physical situation as the interaction of N 
spin-i particles with classical magnetic fields. Denoting 
by Si the spin of ith particle, one can use the collective 
spin operators: 

N 

J = ^S,. (27) 

The orthonormal basis { \j, fi) } consists of the symmetric 
Dicke states Q : 

i^'/^)= j Eni+)'^ni-)'' (28) 
fc=i 

where |-|-); and |— ); are the upper and lower states, re- 
spectively, of the Ith atom, and the summation is over 
all possible permutations of N atoms. If only symmetric 
states are considered, then the "cooperative number" j is 
equal to N/2 and p = ^ + j is just the number of excited 
atoms. Usually the atoms (ions) are far enough apart so 
their wave functions do not overlap and the direct dipole- 
dipole coupling or other direct interactions between the 
atoms may be neglected. 

In the spin formulation (see, e.g., Ref. for a very 
good description), the magnetic moment ^ — iiqS is as- 
sociated with each particle. If a uniform external mag- 
netic field Bo ~ Bqz is applied, the Hamiltonian for each 
particle is given by 

Ha = -fi-Bo = hujoSz, (29) 

where Hujq = ^(IqBq is the separation in energy between 
the two levels. The corresponding Heisenberg equation 
for the collective spin operator is 

d3/dt = LOoX^, (30) 

where u^q = ujqz. Then one applies the so-called clock 
radiation which is a classical field of the form 

'B± = B± {y cosuit — xsinojt) , (31) 

where w 2± and we assume > 0. In the reference 
frame that rotates at frequency u, the collective spin J 
interacts with the effective field 

B = BrZ + B^y, (32) 

where Br — Bo{ujo — uj)/luo. In the rotating frame, the 
Hamiltonian is H — —iiqJ ■ B, and the Heisenberg equa- 
tion for J is 

d3/dt = uj'xJ, (33) 



where u:' = {loq — uj)i -\- uj±y and uj± = —j.ioB±/h. 

The Ramsey method breaks the evolution time of the 
system into three parts. In the first part B± is nonzero 
and constant with value Bi during the time interval 
< t < t^. During this period (the first Ramsey pulse), 
B = BrZ + Biy ~ i?iy, where we made the assumption 
|i?i| 3> \Br\^ i.e., Iwil » |cjo — with wi = -^oBi/h. 
Therefore, in the rotating frame of Eq. (^3|), J rotates 
around the y axis by the angle i?i = uJit^. During the 
second period, of duration T, (usually T :s> t^), B± is 
zero, so B = Br-z, and J rotates around the z axis by 
the angle if — {ujq — uj)T. The third period is exactly as 
the first one but with the field B± — B2 and the corre- 
sponding angular frequency uj2 — —^J■oB2/fl■ This gives a 
rotation around the y axis by the angle 'i?2 = ^2^^. These 
three Ramsey pulses provide the rotations we described 
in Sec. Ill A (usually, i9i = i?2 = 7r/2). 

The aim of spectroscopic experiments is to measure 
the transition frequency luq (which is equivalent to the 
measurement of tp, as uj and T are determined by the ex- 
perimenter). Usually, one measures the number of atoms 
in the upper state |-|-), 

N+out = ^.out + N/2, (34) 

and thus obtains the information about the angle ip or, 
equivalently, about the frequency loq. Of course, in or- 
der to infer this information one should know the initial 
quantum state of the system. The measurement sensi- 
tivity, as seen from Eq. (|25|), also depends on the initial 
quantum state. 

In the state-reconstruction procedure, the first Ram- 
sey pulse should be omitted, while the second and 
third pulses produce the desired phase-space displace- 
ment (?^(n). After the displacement is completed, one 
should measure the probability to find the system in one 
of the states |j, /i), for example, measure the population 
of the ground state \j, —j) or of the most excited state 
\j, j). This measurement can be made by driving a dipole 
transition to an auxiliary atomic level and then observing 
the resonance fluorescence. The phase space is scanned 
by repeating the measurement with many identically pre- 
pared systems for various durations of the Ramsey pulses, 
T and t§. Of course, the apparatus should be first cali- 
brated by measuring the transition frequency luq. 

V. TWO-MODE LIGHT FIELDS 

The basic device employed in a passive optical in- 
terferometer is a beam splitter (a partially transparent 
mirror). A Mach-Zehnder interferometer consists of two 
beam splitters and its operation is as follows. Two light 
modes (with boson annihilation operators oi and 02) are 
mixed by the first beam splitter, accumulate phase shifts 
tpi and tp2, respectively, and then they are once again 
mixed by the second beam splitter. Photons in the out- 
put modes are counted by two photodetectors. In fact. 
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a Michelson interferometer works in the same way, but 
due to its geometric layout the two beam splitters may 
coincide. 

Each beam splitter has two input and two output 
ports. Let a — (ai, 02)"^ and b = (61, 62)"^ be the column- 
vectors of the boson operators of the input and output 
modes, respectively. Then, in the Heisenberg picture, the 
action of the beam splitter is described by the transfor- 
mation 



b = Ba, 



(35) 



where B is a 2 x 2 matrix. For a lossless beam splitter 
B must be unitary, thereby assuring the energy (photon 
number) conservation. A possible form of B is 



cos(i?/2) -sin(i?/2) \ 
sin(??/2) cos(i?/2) ) 



(36) 



with T = cos2(i?/2) and R = sin^(i?/2) being the trans- 
mittance and reflectivity, respectively. When the two 
light modes accumulate phase shifts (fi and ip2, respec- 
tively, the corresponding transformation is 



Pa, 











(37) 



The group-theoretic description of the intcrferometric 
process is based on the Schwinger realization of the 
su(2) algebra: 



= {a\a2 + 4ai)/2, 
Jy = -i{a\a2 - a2ai)/2, 
Jz = {a\ai — a\a2)/2. 



(38) 



Actions of the interferometer elements (mixing by the 
beam splitters and phase shifts) can be represented as 
rotations of the column-vector J = Ux, Jy, Jz)^ ■ The 
beam-splitter transformation of Eq. (|3^) is represented 
by rotation Ry(t?) around the y axis by the angle "d, 
and the phase shift of Eq. ( |37| ) is represented by rota- 
tion Rz(<p) around the z axis by the angle ip = ip2 — (pi- 
Now, if the transmittances of the two beam splitters are 
Ti = cos2(i9i/2) and T2 = cos'^{-d2/2), respectively, then 
the interferometer action is given by the three rotations 
described in Sec. Ill A. (Usually, one uses 50-50 beam 
splitters, so i9i = —■d2 — 7r/2.) 

Interferometers are constructed to measure the relative 
phase shift (/p, which is proportional to the optical path 
difference between the two arms. Usually, one measures 
the difference between the photocurrents due to the two 
output light beams. This quantity is proportional to the 
photon-number difference at the output, gout = 2 J^out- If 
the input state of light is known, then the measurement 
of Qout can be used to infer the phase shift ip and estimate 
the measurement error due to the quantum fluctuations 
of the light field. 

A simple calculation gives — (iV/2)(l + A^/2), where 
N = a\ai + a\a2 is the total number of photons in the 



two modes. If N has a fixed value for the input state 
of the two-mode light field, then this state belongs to 
the Hilbert space Tij of a specific SU(2) representation 
with j = N/2. Because N is the SU(2) invariant, this 
state will remain in Tij during the intcrferometric pro- 
cess. Such input states of the two-mode light field can be 
reconstructed using a rearrangement of the intcrferomet- 
ric scheme, according to the general procedure described 
in Sees. II and III B. 

The phase-space displacement (7^(n) needed for the 
state-reconstruction procedure can be implemented by 
using an interferometer without the first beam splitter. 
Then one should measure the probability P;^(n) to find 
the output light in one of the states |j, A*)- Note that 
these states are given by 



|j» = b + m)i ® li - M)2 



(39) 



in the terms of the Pock states of the two light modes. 
So, /i is just one half of the photon-number difference 
measured at the output. Averaging over many measure- 
ments, one obtains the probabilities p^i{'n)- For example, 
Pj(n) is the probability that all photons exit in the first 
output beam while the number of photons in the sec- 
ond output beam is zero. The measurement should be 
repeated with identically prepared input light beams for 
many phase-space displacements. This means that one 
needs a well-calibrated apparatus which can be tuned for 
various values of the relative phase shift ip. These phase 
shifts can be conveniently produced by moving a mirror 
with a precise electro-mechanical system. Various values 
of the angle •§2 can be realized using a collection of par- 
tially transparent mirrors with different reflectivities for 
the second beam splitter. An alternative possibility is 
to use the dependence of the reflectivity on the angle of 
incidence for light polarized in the plane of incidence. 

In general, the state reconstruction for two-mode light 
fields is a tedious task, because the corresponding Hilbert 
space is very large pl|-p5[. Obviously, this task can be 
greatly simplified for the subclass of two-mode states 
with a fixed total number of photons, by means of the 
reconstruction method presented here. However, this 
method is in principle suitable also for other two-mode 
states as well. In general, the whole Hilbert space of the 
two-mode system can be decomposed as 



(40) 



The method of inverted interferometry enables one to re- 
construct the part of the density matrix corresponding 
to each irreducible subspace TLj . One case for which our 
method is applicable is the subclass of states, whose den- 
sity matrices are block-diagonal in terms of the decompo- 
sition (|^). This means that the corresponding operator 
can be written as 



(41) 
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where pj is an operator on Tij. Each component pj 
evolves independently during the phase-space displace- 
ment; hence the state of the whole system can be mea- 
sured by reconstructing all invariant components pj . The 
other case for which our method works is the subclass of 
pure states, 

3 

\^) = E I^J-)' = E c.mIj'.m)- (42) 

Then the density matrix can be written as 

p = EiV'.-)(^^i + Ei^^)(^^'i- (43) 

3 3 ¥=3' 

The populations of the states \j, p) are unaffected by the 
second term in (^ ) , and one can reconstruct all invariant 
components pj = \tpj){ipj\. This gives information about 
the state of the whole system, except for relative 
phases between different From the technical point of 
view, each measurement of the photon-number difference 
2p, needed to determine the probabilities P/j(n), should 
be accompanied by a measurement of the photon- number 
sum N = 2j, in order to determine to which invariant 
subspace Tij does the detected value of p correspond. 
Consequently, one needs to make many more measure- 
ments, in order to accumulate enough data for each value 
of J. A technical problem is that quantum efficiencies of 
realistic photodetectors are always less then unity. While 
this problem is not too serious for the measurement of 
the photon-number difference (as long as both detectors 
have the same efficiency), it puts a serious limitation on 
the accuracy of the measurement of the total number of 
photons. 



VI. TWO-DIMENSIONAL VIBRATIONS OF A 
TRAPPED ION 

As was recently demonstrated by Wineland et al. , 
a single laser-cooled ion in a harmonic trap can be used to 
simulate various interactions governing many well-known 
optical processes. In particular, one can simulate trans- 
formations produced by elements of a Mach-Zehnder op- 
tical interferometer. 

Consider a single ion confined in a two-dimensional 
harmonic trap, with angular frequencies of oscillations 
in two orthogonal directions f2i and Q2- Two internal 
states of the ion, |+} and |— ), are separated in energy 
by hujQ. The internal and motional degrees of freedom 
can be coupled by applying classical laser beams, with 
electric fields of the form 

E(x, t) = Eo cos(k ■ X - ojt + 

For example, one can apply two laser beams to pro- 
duce stimulated Raman transitions. We denote by a; = 



LU1—UJ2, k = ki — k2, and <I> = $1 — $2 the differences be- 
tween the angular frequencies, the wave vectors, and the 
phases, respectively, of the two applied fields. Then, in 
the rotating-wave approximation, the interaction Hamil- 
tonian reads 

i// = nKexp[i(k-x-(5t + $)] +H.C., (44) 

where S = lu — luq is the frequency detuning, x is the ion's 
position relative to its equilibrium, and k is the coupling 
constant (the Rabi frequency). Each of the two modes of 
the ion's motion can be modelled by a quantum harmonic 
oscillator: 

X,. = XOr{ar + al), XQr ~ y/hJJ^MTlr) , (45) 

where r = 1,2 and M is the ion's mass. Also, let 
rjr — krXor (r — 1,2) be the Lamb-Dicke parameters 
for the two oscillatory modes. It is convenient to use the 
interaction picture for the ion's motion: 

Hj = exp{iHot/h)Hi exp{—iHot/h) 

= fce*^*-'^*) Y[ exp[ir]ridr + 4)] + H-C, (46) 

r=l,2 

where Hq is the free Hamiltonian for the ion's motion. 

Ho = Mil {a\ai + 5] + fif^2 (4^2 + 5) , (47) 

and CLr = ttr exp(—iflrt) , r = 1,2. 

If the coupling constant k is small enough and 
and Q2 are incommensurate, one can resonantly excite 
only one spectral component of the possible transitions. 
For a particular resonance condition 6 — — ^1 (and 
in the Lamb-Dicke limit of small 771 and 772), the prod- 
uct in Eq. ( ^6|) will be dominated by the single term 
{iTiiai){iri2a2}- Therefore, one obtains 

Hi w -hmnm (e'*ai4 + 6-^*402) . (48) 

Returning to the Schrodinger picture, the total evolution 
operator reads: 

U{t) = cx-p{-iHot/h) cxp{-iHit/h) 

= exp[-z(rii + n2)iN + l)t/2] cxp[i{n2 - ^i)Jzt] 
X exp(2zK77i772 J$<). (49) 

Here, N — a\ai +a2a2 is the total number of vibrational 
quanta in the two modes, J$ = Jx cos $ -f Jj, sin <I>, and 
we used the Schwinger realization ( ^8|) for the SU(2) gen- 
erators. 

Now, let us consider only such motional states of the 
ion for which N has a fixed value, i.e., which belong to the 
irreducible Hilbert space Hj (with j = N/2). For these 
states, the first exponent in (|4^ ) will just produce an 
unimportant phase factor and can be omitted. Clearly, 
the evolution operator ( ^9| ) can be used to simulate the 
action of an optical interferometer, with two vibrational 
modes of a trapped ion employed instead of two light 
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beams. In order to simulate the action of a beam split- 
ter, one should apply the interaction (^) during time tg 
and ensure that \2Kr]ir]2\ » |ri2 — so the effect of the 
free evolution can be neglected. Then, for <i> = 7r/2, the 
evolution operator reads 



Uy{9) = exp{i6Jy), 9 — 2Kr]ir]2t0. 



(50) 



A relative phase shift between the two modes can be 
produced just by using the free evolution, i.e., with no 
external laser fields applied. Letting the system evolve 
freely during time T, one obtains 

C/,(0) =exp(z0J,), ^={n2-ni)T. (51) 

It is obvious that applying consequently the transfor- 
mations ( |5l|) and ( ^o| ) one will produce the phase-space 
displacement 5^(n), employed in the state-reconstruction 
procedure. The whole phase space can be scanned by re- 
peating the procedure with identically prepared systems 
for various durations T and tg. Each phase-space dis- 
placement should be followed by the measurement of the 
probability (n) to find the system in one of the states 
\j, j-i). For example, Pj(n) is the probability that the first 
oscillatory mode is excited to the iVth level {N = 2j) 
while the second mode is in the ground state. Such a mea- 
surement can be made with the method used recently by 
the NIST group ||2^] to reconstruct the one-dimensional 
motional state of a trapped ion. The principle of this 
method is as follows. One of the oscillatory modes is cou- 
pled to the internal transition |-|-) ^ |— ). This is done by 
applying one classical laser field, so single-photon tran- 
sitions are excited. This results in an interaction of the 



Jaynes-Cummings type |28| between the oscillatory mode 
and the internal transition. Then the population P_ (i) of 
the lower internal state |— ) is measured for various values 
of the interaction time t (as we already mentioned, this 
measurement can be made by monitoring the resonant 
fluorescence produced in an auxiliary dipole transition). 
If |— ) is the internal state at < = 0, then the signal aver- 
aged over many measurements is 



n=0 



F„ cos(2fi„_„+it)e 



where fi„_„_(_i are the Rabi frequencies and 7„ are the ex- 
perimentally determined decay constants. This relation 
allows one to determine the populations P„ of the mo- 
tional eigenstates \n). By virtue of Eq. (|39|), this gives the 
populations of the SU(2) states |j, ^) (with fi = n — j 
for the first mode and fi — j — n for the second mode). 
For example, p-j and pj are given by Pq for the first and 
second modes, respectively. 



VII. CONCLUSIONS 

In this paper we presented practical methods for the re- 
construction of quantum states for a number of physical 



systems with SU(2) symmetry. All these methods em- 
ploy the same basic idea — the measurement of displaced 
projectors — which in principle is applicable to any sys- 
tem possessing a Lie-group symmetry. Practical realiza- 
tions, of course, vary for different physical systems. In 
our approach, we exploited the fact that transformations 
applied in conventional spectroscopic and interferomet- 
ric schemes are, from the mathematical point of view, 
just rotations. In the context of the SU(2) group, these 
rotations constitute phase-space displacements needed to 
implement a part of the reconstruction procedure. There- 
fore, the spectroscopic and interferometric measurements 
can be easily rearranged in order to enable one to deter- 
mine unknown quantum states for an ensemble of iden- 
tically prepared systems. As the spectroscopic and in- 
terferometric measurements are known for their high ac- 
curacy, we hope that the corresponding rearrangements 
will allow accurate reconstructions of unknown quantum 
states. 
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